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Resumo 

O saldo de radiação em ecossistemas possui comportamento complexo e representa a energia disponível para os processos físicos e bioló- 
gicos que ocorrem entre a superfície vegetada e a atmosfera. Sendo assim, o objetivo deste trabalho foi analisar a complexidade do saldo 
de radiação medido acima do dossel de uma floresta sazonalmente inundada no Norte do Pantanal, sob a ótica da teoria da complexi- 
dade. Para tanto, foram reconstruídos os atratores sazonalmente, estimando os expoentes dimensão de correlação (E> 2(chuvoso) = 1,73±0,05 e 
D 2 ( sec o) = 1,36±0,18), dimensão embutida (m = 7), taxa de recorrência (RR (dmvoso) = 0,144±0,004 e RR (seco) = 0,169±0,002) e o determinismo 
(DET, , , = 0,896±0,013 e DET, , = 0,973±0,015) das séries de dados analisados. Os resultados evidenciam a existência de atratores 

estranhos de baixa dimensão fractal e sensíveis à sazonalidade do saldo de radiação. 

Palavras-chave: atratores, gráficos de recorrência, dinâmica não linear. 


Abstract 

The ecosystem net radiation has a complex variation, and represents the energy available for physical and biological processes that occur 
between the vegetated surface and atmosphere. Therefore, the objective of this study was to analyze the complexity of net radiation measu- 
red above the canopy of a seasonal flooded forest in the northern Pantanal from the perspective of complexity theory. The attractors were 
seasonal reconstructed by the estimation of the exponents correlation dimension (D 2(wet) = 1,73±0,05 and E) 2(dry) = 1,36±0,18), embedded 
dimension (m = 7), rate of recurrence (RR (wet) = 0,144±0,004 and RR (dry) = 0,169±0,002) and determinism (DET (wet) = 0,896±0,013 and 
DET (dry) = 0,973±0,015) data sets analyzed. The results showed the existence of strange attractors of low fractal dimension, and sensitive 
to the seasonality of radiation balance. 

Keywords: attractors; recurrence plot; nonlinear dynamics. 
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INTRODUÇÃO 

A complexidade dos ecossistemas flores- 
tais (tropicais, inundados e/ou cerrado), dado o 
seu número e a variedade de seus subsistemas em 
interação, assim como os efeitos antropogênicos, 
torna a sua descrição um problema brasileiro extre- 
mamente significativo. (Gomes & Varriale, 2004). 

O Pantanal é a maior planície tropical de 
inundação sazonal do planeta, com aproximada- 
mente 138 mil km 2 , e abriga biodiversidade despro- 
porcionalmente alta e influencia substancialmente 
o ciclo de carbono e outros gases do efeito estufa. 
O Pantanal distingue-se de outros biomas por 
apresentar fauna e flora peculiar que desempenha 
papel muito importante no saldo de radiação, na 
manutenção dos fluxos de energia e vapor de água, 
regulação da temperatura e estabilidade dinâmica 
do clima regional e global (Silva & Abdon, 1998; 
Millenium Ecosystem Assesssment, 2005; Junk & 
Nunes da Cunha, 2005). 

Pesquisadores afirmam que alterações no 
ciclo hidrológico, observado pelo aumento dos 
picos de precipitação máxima, sugerem que a 
região está sendo afetada pelas mudanças climáticas 
(Diniz Fonseca, M.; Campeio Jr et al., 2008). Com 
baixa drenagem, o solo pantaneiro frequentemente 
prolonga o período de alagamento, limitando a 
estação seca de abril até setembro. A seca causa 
estresse hídrico nas plantas locais, regulando sua 
fenologia e produção de fitomassa (Junk, 2002). Há 
ainda rápida mudança na vegetação em campos 
sazonalmente inundados que também tem sido 
considerado um efeito de mudanças climáticas 
ocorridas em ciclos plurianuais (Arieira & Nunes 
da Cunha, 2006). O Pantanal é suscetível a uma 
série de consequências ecológicas causadas pelo 
desmatamento, como o aumento da temperatura do 
ar, redução da precipitação e o prolongamento da 
seca pela ampliação das áreas destinadas à criação 
de gado e agricultura na região (Nobre et al., 1996; 
McClain, 2002; Costa & Pires, 2010; Lewis et al., 
2011, Lapola et al, 2011). 

Com a troca da cobertura vegetal original 
por áreas de pastagens e cultivos, o fluxo de vapor 
de água para a atmosfera é alterado, bem como a 
refletividade do solo, afetando diretamente o saldo 
de radiação entre outros efeitos. Biomas tropicais 
são sistemas abertos, dinâmicos e complexos, o 
que os torna importantes para a avaliação de seus 
balanço de energia, no que diz respeito às estima- 
tivas da evapotranspiração (ET), aquecimento do 
ar e do solo, e ainda prever as suas contribuições 
com as mudanças climáticas locais e globais a 
médio e longo prazo (Pereira et al., 2002; IPCC, 
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2007, Motta et al., 2011). 

Ainda não se compreende como essas alte- 
rações regionais no uso do solo podem afetar o 
ciclo hidrológico e o regime de chuvosos da região. 
Um dos motivos é a falta de séries históricas de 
dados micrometeorológicos medido in situ para 
detectar e monitorar estas possíveis mudanças no 
clima regional. 

Embora haja estudos das relações entre 
radiação líquida e suas partições no Pantanal, 
ecossistemas são sistemas não determinísticos, 
não lineares, em que a influência de um estado 
anterior da dinâmica sobre o estado posterior não 
é diretamente proporcional, mas envolve laços de 
retroalimentação que podem influenciar o sistema 
globalmente (Souza & Buckeridge, 2004). Além 
disso, a recorrência de estados, no sentido de que 
os estados arbitrários do sistema poderão estar 
próximos novamente com o avanço do tempo, é 
uma propriedade fundamental de sistemas dinâ- 
micos determinísticos e é típico para sistemas não 
lineares (Marwan, 2007). 

Sistemas abertos fora do equilíbrio são 
caracterizados por atratores, que são regiões restritas 
do espaço de estados para onde as trajetórias do 
sistema convergem. No caso de um sistema com 
dinâmica puramente aleatória, o atrator do sistema 
preenche o espaço de estados inteiramente, isto 
é, não há nenhum tipo de organização temporal, 
a dinâmica do sistema é totalmente imprevisível 
a curto ou longo prazo. Enquanto que em um 
sistema aberto e fora do equilíbrio sua dinâmica 
fica restrita a uma região do espaço de estados do 
atrator (Souza & Buckeridge, 2004). 

A importância da aplicação de ferramentas 
adequadas para se compreender os processos de 
troca de energia, matéria e momentum em floretas 
tropicais é destacada em pesquisas que abordam 
a natureza não linear dos dados ambientais, con- 
tribuindo para o entendimento e modelagem da 
interação solo-planta-atmosfera. 

Embasado na teoria dos sistemas dinâmicos, 
várias pesquisas envolvendo séries temporais de dife- 
rentes variáveis ambientais mencionam a existência 
de atratores climáticos de baixa dimensão fractal e 
indicam o número mínimo de variáveis presentes 
na evolução da dinâmica do sistema (Ruelle & 
Takens,1971; Nicolis & Nicolis, 1984; Fraedrich, 
1986; Lorenz, 1991; Waldrop, 1992; Povedo-Jara- 
millo & Puente, 1993; Nicolis & Prigogine, 1998; 
Xin et al., 2001; Gallego et al., 2001; Peitgen et al., 
2004; Campanharo et al., 2008). 

Por definição, um Gráfico de Recorrência 
(GR) é uma técnica avançada de análise de dados 
de sistemas não lineares por meio da visualização 
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de gráficos formados por uma matriz quadrada, 
em que os elementos da matriz correspondem aos 
momentos em que um estado do sistema dinâmico 
se repete (colunas e linhas correspondem então a 
um determinado par de vezes). Tecnicamente, o 
GR revela todos os momentos em que a trajetória 
do espaço de fase do sistema dinâmico visita mais 
ou menos a mesma região no espaço de fase. 

O estudo dos GRs dos sistemas dá uma 
boa introdução sobre a tipologia característica e 
a textura. No entanto, sua quantificação oferece 
uma maneira mais objetiva para a investigação do 
sistema considerado. As técnicas de quantificação 
para análise de dados, ou seja, a quantificação das 
medidas de complexidade dos GRs, têm se tornado 
muito popular num grupo crescente de cientistas 
de diversas áreas. Sendo assim, o objetivo deste 
trabalho foi analisar a complexidade do saldo de 
radiação medido acima do dossel de uma floresta 
sazonalmente inundada no Norte do Pantanal, 
sob a ótica da teoria da complexidade. Para tanto, 
foram utilizados GR através da representação 
bidimensional de trajetórias no espaço de fase (o 
que é vantajoso por se tratar de um sistema de alta 
dimensão), sendo possível a verificação dos padrões 
típicos que estão ligados ao comportamento espe- 
cífico do sistema; e as Medidas de Quantificação 
de Recorrência (MQR) para contabilização dos 
pontos e das diagonais, quantificando a complexi- 
dade pela taxa de recorrência (RR) e determinismo 
(DET) do sistema. 

Gráficos de Recorrência 

O gráfico de recorrência (GR) representa o 
número de vezes que o sistema recorre a um estado 
anterior, ou seja, em que uma trajetória visita a 
mesma região no espaço de fase. Essa trajetória no 
espaço de fase pode ser reconstruída a partir de 
uma única série temporal pela incorporação de 
tempos de atrasos (x), segundo o teorema de Talcens 
(1981). O GR é obtido após algumas etapas. Pri- 
meiramente, obtém-se (ou reconstrói-se) outras 
variáveis (X k (t) } (sendo k=l, ..., n-1) conhecendo 
somente X Q (t). Depois, desmembra-se a série tem- 
poral original X Q (t) em deslocamentos temporais 
sucessivos de defasagem fixa x (x = mAt, onde m é 
um número inteiro e At é um intervalo entre suces- 
sivas amostras) para “n” pontos equidistantes do 
conjunto de dados. Isto é: 

Xo: Xo(ti), . . . , Xo(Je) 

Xi! Xo(ti + x), . . . , Xo(ÍN + t) 

X n - 1 : X 0 [ti + (n - 1 ) x] 5 . . . , X 0 [tN + (n - 1) x] 


em que a melhor escolha de x deve apresentar uma 
independência linear entre a série original e as séries 
defasadas, ou seja, a menor correlação, no intuito 
de obter a maior informação possível do sistema 
estudado. 

Eckmann et al. (1987) introduziram uma 
ferramenta capaz de visualizar a recorrência de 
estados x em um espaço de fase. Geralmente, um 
espaço de fase tem mais de três dimensões, o que 
não permite que ele seja retratado. Espaços de 
fase com dimensões mais elevadas (hiperespaço) 
só podem ser visualizados pela projeção em sub- 
-espaços de duas ou três dimensões. No entanto, 
a ferramenta de Eckmann et al. (1987) permite 
investigar a trajetória ra-dimensional do espaço de 
fase através de uma representação bidimensional de 
suas repetições. A recorrência de um estado i em um 
tempo diferente j é marcada dentro de uma matriz 
bidimensional quadrada com pontos 1 e 0 (pontos 
pretos e brancos respectivamente), onde ambos os 
eixos são eixos de tempo. Esta representação é o 
próprio Gráfico de Recorrência (GR). Um GR pode 
ser matematicamente expresso como: 

Ri j = 6(s— 1| x t — Xj ||), x t e fi m , i,j = 1 ...N, 

em que N é o número de estados considerados x p 
8 é a distância limite, II * II é a norma e 8 ( * ) é a 
fúnção de Heaviside. Para estados recor- rentes, 
isto é, para estados que estejam em uma região de 
vizinhos em um raio , se introduz a seguinte nota- 
ção (Marwan et al., 2007): x t « Xj <=> R t j = 1 

O Gráfico de Recorrência é obtido pela plo- 
tagem da matriz de recorrência e traçando um ponto 
para as coordenadas (ij) se R = 1 , e um ponto branco 
se R = 0. O GR sempre tem uma linha diagonal 
principal, a linha de identidade. Além disso, o GR 
é simétrico no que diz respeito a diagonal principal. 
A escolha do melhor valor para raio de sondagem 
tem sido sugerido e defendido na literatura como 
uma pequena porcentagem do diâmentro máximo 
do espaço de fase (Mindlin & Gilmore, 1992). 

O propósito inicial de GR é a inspeção visual 
de trajetórias em espaços de fase de altas dimensões. 
A visão de um GR dá sugestões sobre a evolução 
temporal destas trajetórias. A vantagem de GR é 
que eles também podem ser aplicados a conjuntos 
de dados curtos e não estacionários. 

Os gráficos de recorrência exibem carac- 
terísticas de grande escala e padrões de pequena 
escala. Os primeiros padrões de forma denotados 
por Eckmann et al. (1987) foram a tipologia e a 
textura. 
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Figura 1 - (A) Segmento da trajetória no espaço de fase do sistema de Lorenz (parâmetros de controle r = 28, a = 10, b = 03/08; 
Lorenz, 1963), usando as três componentes x, y e z; e (B) correspondente a sua parcela de recorrência. Um ponto da trajetória de 
j dentro da vizinhança (círculo cinza em (A)) de um determinado ponto no i é considerado como um ponto de recorrência (ponto 
preto sobre a trajetória de (A)). Este é marcado com um ponto preto no RP no local (i, j). Um ponto fora da vizinhança (pequeno 
círculo de (A)) provoca um ponto branco no RP. O raio da vizinhança para o GR é e = 5 (Disponível em < http://www.recurrence- 
-plot.tk/glance.php > acesso em 12/2011). 



Figura 2 - Gráficos de Recorrência de (A) movimento periódico de uma frequência, (B) atrator estranho de Rossler (parâmetros 
a=b=0.2 e c=5.7) e (C) ruído branco (adaptado de Marwan et al., 2007). 


A tipologia oferece uma impressão global 
que pode ser caracterizada como homogênea, flu- 
tuante, periódica e desviante. Sistemas homogêneos, 
típicos de sistemas estacionários e autônomos, 
possuem tempo de relaxação curto em comparação 
com o tempo abrangido pelo GR. Um exemplo é 
o GR de uma série temporal de dados aleatórios. 
Em sistemas oscilantes o GR contêm diagonais 


orientadas, estruturas periódicas recorrentes (linhas 
diagonais, estruturas quadriculadas). 

Em sistemas quase periódicos, as distâncias 
entre as linhas diagonais no GR são diferentes. 
No entanto, mesmo em sistemas oscilantes cujas 
oscilações não são facilmente reconhecíveis, os 
GRs podem ser utilizados a fim de encontrar 
suas oscilações. Sistemas não estacionários são 
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os desvios causados por sistemas com parâmetros 
de variação lenta. Tal mudança lenta (adiabática) 
ilumina o canto superior esquerdo e inferior direito 
do GR. Mudanças abruptas na dinâmica, bem 
como eventos extremos causam áreas brancas ou 
bandas no GR. 

Uma inspeção mais refinada dos GRs revela 
estruturas de pequena escala (a textura), que são 
pontos isolados, linhas diagonais, bem como linhas 
verticais e horizontais (a combinação de linhas 
verticais e horizontais, formas retangulares aglo- 
merados de pontos de recorrência). Podem ocorrer 
pontos de recorrência isolados se os estados são 
raros, se não persistir por muito tempo ou se eles 
oscilam muito. No entanto, eles não são um sinal 
exclusivo do acaso ou de ruído (por exemplo, em 
mapas de atratores estranhos). 

A linha diagonal R. +k j+k - 1 (para k = 1 ... /, 
onde / é o comprimento da linha diagonal) ocorre 
quando um segmento da trajetória segue paralelo a 
outro segmento, ou seja, a trajetória visita a mesma 
região de espaço a fase em momentos diferentes. 
O comprimento desta linha diagonal principal é 
determinado pelo tempo de duração semelhante 
da evolução local dos segmentos de trajetória. 

A direção dessas estruturas diagonais pode 
ser diferente. Linhas diagonais paralelas à linha 
principal representam trajetórias paralelas seguindo 
uma mesma evolução temporal. As estruturas 
diagonais perpendiculares à linha diagonal prin- 
cipal representam o um seguimento de evolução 
contrária. 

A linha vertical (horizontal), Ri j + k = 1 
(para k = 1 ... v, onde v é o comprimento da linha 
vertical) marca um período de tempo em que um 
estado não muda ou muda muito lentamente. Ao 
que parece, estes estados são presos por algum 
tempo, comportamento típico de estados laminares 
(intermitência). Essas estruturas de pequena escala 
são a base de uma análise quantitativa dos GRs. 

Análise da Complexidade 

Embora as estruturas GR possam por si 
só caracterizar alguns sistemas, dando forma às 
diferentes tipologias, elas desempenhas um papel 
fundamental nas quantidades relacionadas ao 
padrão de pequena escala, a chamada textura. 
Estas quantidades se baseiam tanto na densidade 
de recorrência quanto nas quantidades e tamanhos 
das linhas diagonais e verticais (horizontais). 

A análise de complexidade é dada pelas 
Medidas de Quantificação de Recorrência (MQR), 
que é um método de análise de dados não lineares 
que quantifica o número e duração das recorrên- 


cias de um sistema dinâmico representado por 
sua trajetória no espaço de fase das variáveis. A 
análise da MQR foi desenvolvida por Zbilut & 
Webber Jr (Zbilut & Webber Jr., 1992; Webber Jr 
& Zbilut, 1994) e estendida com novas medidas 
de complexidade por Marwan et al. (2002). 

MATERIAL E MÉTODOS 

Sítio Experimental e Instrumentação 

A área experimental localiza-se na Reserva 
Particular do Patrimônio Natural - RPPN SESC - 
Pantanal, município de Barão de Melgaço - MT, 
distante 160 km de Cuiabá, Mato Grosso, na Região 
Centro-Oeste do Brasil onde está instalada uma torre 
micrometeorológica de 32 m de altura (16°39’50”S; 
56°47’50”O, altitude de 120 m). 

A torre é rodeada por uma vegetação flores- 
tal monodominante de cambará ( Vochysia divergem), 
localmente chamada de cambarazal, em uma faixa 
contínua de 25 km de extensão por 4 km de largura, 
paralela ao rio Cuiabá (Junk, 2002) e com altura do 
dossel entre 28 a 30 m. O cambarazal tem apresen- 
tado rápido espalhamento nos campos inundados 
sazonalmente (Arieira & Nunes da Cunha, 2006). 

O solo é GLEISSOLO HÁPLICO Ta 
Distrófico. O clima é AW segundo Kõppen, a 
temperatura do ar média anual oscila entre 22°C 
e 32°C e a precipitação média anual entre 1 100 e 
1200 mm com período seco de maio a setembro 
(Hasenack, 2003). 


Tabela 1 - Subdivisão dos dados por mês/ano e características 
sazonais (período chuvoso e seco) em uma floresta sazonalmen- 
te inundada no norte do Pantanal 


Mês / Ano 

Sazonalidade 

Setembro de 2006 

Seco 2006 

Fevereiro de 2007 

Chuvoso 2006/07 

Agosto de 2007 
Janeiro de 2008 

Seco 2007 
Chuvoso 2007/08 

Agosto de 2008 
Janeiro de 2009 

Seco 2008 

Chuvoso 2008/09 

Agosto de 2009 
Nov/Dez de 2009 

Seco 2009 
Chuvoso2009/10 


Os valores de saldo de radiação, analisados 
neste trabalho, foram medidos por meio de um saldo 
radiômetro (Net Radiometer, Kipp & Zonen Delft, 
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Inc., Holland) instalado a 33 m de altura. Os dados 
foram coletados a cada 1 segundo e armazenados 
em médias de 15 minutos em um datalogger (CR 
10X, Campbell Scientific, Inc., Ogden, Utah). 

Os dados desse trabalho são referentes ao 
período seco de 2006 até o período chuvoso de 
2009/2010, perfazendo um espectro de oito (8) 
sazonalidades, e aproximadamente quatro (4) anos 
(Tabela 1). 


Métodos 

Neste trabalho, foram utilizadas duas Medi- 
das de Quantificação de Recorrência, a taxa de 
recorrência (RR) e o determinismo (DET). 

A taxa de recorrência (RR) mede a densi- 
dade de pontos de recorrência no GR, definido 
pela equação abaixo. 


N 

^ X Rl > (£) 

i,j = i 

A RR não exclui a diagonal principal e, no 
limite N~* °o é a probabilidade de que um estado 
retorne à sua vizinhança de raio 8 no espaço de fase. 

O determinismo {DET) é a medida de 
previsibilidade do sistema, definido pela equa- 
ção abaixo. 


DET = 


Zi-i ■ IP ( 0 

SU ip(1) 


em que P(l) é o número de diagonais de com- 
primento/ no GR, e / é o menor tamanho 
para uma linha ser considerada uma diagonal. 

Os cálculos foram realizados no Matlab®, 
utilizando um pacote definido como TOCS Y - Tool- 
box for Complex Systems pelo Potsdam Institute 
for Climate Impact Research (PIK) (disponível 
em < http://tocsy.agnld.uni-potsdam.de > acesso 
em 27/04/2011). 


RESULTADOS E DISCUSSÃO 

A reconstrução do atrator com coordenadas 
de atraso temporal não é um processo automático, 
logo é importante a escolha de um x apropriado 
(Baker & Gollub, 1996). Palú (2008) usou o método 
da Informação Mútua (Abarbanel et al., 1993; Baker 


& Gollub, 1996) para escolher o melhor tempo de 
defasagem temporal (x) em análise de dados prove- 
nientes da torre micrometeorológica na Amazônia 
Legal em Sinop - MT. Palú (2008) afirma que um 
tempo de defasagem x = 5 horas pode ser adotado em 
quaisquer que sejam as análises de séries temporais 
e para a reconstrução de atratores dos processos 
envolvidos na interação atmosfera-biosfera daquela 
floresta, independentemente de qual variável estiver 
sendo estudada e de que época do ano os dados se 
referem não existindo diferenciação entre as estações. 
Assim, definiu-se para a variável uma defasagem 
temporal (x) igual a 5 h (mAt = 20x15 min, em que 
m é um número inteiro e At é um intervalo entre 
cada amostra), visto a baixa correlação entre a série 
original e as série defasada. Os resultados de Palú 
(2008) corroboram os de Baker & Gollub (1996) em 
que a escolha do parâmetro x, em situações práticas, 
pode ser regulado até que os resultados se tornem 
satisfatórios. Também que uma pequena variação 
no valor de x não é significativa para alterar os 
resultados obtidos pelo método de análise. 

Neste trabalho, também foi calculado o 
atraso x pela função auto-correlação como pro- 
posto por Kantz & Schreiber (1996), encontrando 
o mesmo valor de 5 h. O melhor espaço de imersão 
m para evitar a interferência dos falsos vizinhos foi 
calculado pelo algoritmo de Grassberger & Pro- 
caccia (1983), assim como o expoente dimensão 
de correlação (D2), valor aproximadamente igual 
à dimensão fractal do atrator do sistema. 

A estimativa para a dimensão de correlação 
do balanço de radiação pelo algoritmo de Grass- 
bergue & Procaccia (1983) foi de E) 2(chuvoso) = 1,73 
± 0,05 e E) 2(seco) = 1,36 ± 0,18) e para a dimensão 
embutida m - 7 para as duas sazonalidades (Figura 
3 e Figura 4). 

Existe um patamar de saturação alcançado 
em quase todas as análises e observa-se em algumas 
um pequeno acréscimo na dimensão de correlação, 
o que pode ser devido a alguma componente não 
linear no sistema. 

Para melhor visualizar e analisar as trajetó- 
rias no espaço de fase, verificou-se os seus padrões 
típicos e as suas ligações com o comportamento 
específico do sistema, todos os GRs (Figura 5 e 
Figura 6) foram calculados reconstruindo atratores 
embutidos em um espaço de fase de 7 dimensões 
(m = 7), com um atraso (x) de 5 horas, usando a 
norma euclidiana (esfera de sondagem) com s- 2, 
conforme e Mindlin & Gilmore (1992). 

Houve maior flutuação e maiores médias 
do balanço de radiação durante o período chu- 
voso, ocorrendo o oposto durante o período seco. 
Durante os meses de setembro de 2006, agosto de 
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Figura 3 - Estimativa da melhor dimensão de imersão m para o atrator a partir do comportamento da dimensão embutida ra(Rn) 
versus dimensão de correlação D 2 referentes aos períodos secos: (a) 2006, (b) 2007, (c) 2008 e (d) 2009, em uma floresta sazonalmente 
inundada no norte do Pantanal. 
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Figura 4 - Estimativa da melhor dimensão de imersão m para o atrator a partir do comportamento da dimensão embutida ra(Rn) 
versus dimensão de correlação D 2 referentes aos períodos chuvosos: (a) 2006/2007, (b) 2007/2008, (c) 2008/2009 e (d) 2009/2010, 
em uma floresta sazonalmente inundada no norte do Pantanal. 
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(d) 


Figura 5 - Períodos secos. Gráficos de recorrência do balanço de radiação (Rn) referentes à (a) setembro de 2006, (b) agosto de 2007, 
(c) agosto de 2008 e (d) agosto de 2009, em uma floresta sazonalmente inundada no norte do Pantanal. 



Figura 6 - Períodos chuvosos. Gráficos de recorrência do balanço de radiação (Rn) referentes à (a) fevereiro de 2007, (b) janeiro de 
2008, (c) janeiro de 2009 e (d) novembro e dezembro de 2009, em uma floresta sazonalmente inundada no norte do Pantanal. 
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2007, agosto de 2008 e agosto de 2009 (Figura 5) 
apresentaram minúcias de conjuntos de estados 
não estacionários (raros ou muito longe do normal), 
podendo ainda ter ocorrido uma transição de estado 
(Marwan et al.; 2007). Percebe-se comportamento 
cíclico pela distância temporal entre os padrões 
das linhas paralelas que correspondem ao período. 

A evolução dos estados é semelhante em 
diferentes épocas, tendendo a um processo deter- 
minista, visto as linhas diagonais paralelas à linha 
identidade. Breves faixas brancas observadas nas 
figuras 4a, 4b e 4c e levemente em 4d também são 
características de dados não estacionários, ou seja, 
podem ter ocorridos alguns estados diferenciados, 
transições de fase no mecanismo de balanço de 
radiação. Formações de aglomerados ( clusters ) em 
todos os gráficos dos períodos secos são indicado- 
res de estados que não mudam, ou mudam muito 
lentamente no decorrer do tempo, os chamados 
estados laminares (Zbilut & Webber Jr., 1992). 
Esses aglomerados, provavelmente, são relativos 
à presença marcante de massas de ar frio (frentes 
frias) que chegam ao Pantanal pelo sul do Brasil. 
Essas frentes frias, conhecidas como friagens, são 
enormes massas de ar frio que provocam nevoeiro 
e costumam ter a duração de 3 a 7 dias, passando a 
influenciar significativamente a dinâmica atmosfé- 
rica da região (Webber Jr & Zbilut, 1994). Mesmo 
com esses destaques, o padrão dominante dos 
gráficos de recorrência para os períodos secos é 
semelhante ao padrão teórico da figura 2a, movi- 
mento periódico de uma frequência, mas permeado 
de intermitências. Houve uma queimada de elevada 
extensão em 2006 e a fumaça soma-se as causas do 
padrão do gráfico de recorrência. 

Ressalta-se que setembro de 2006 apresenta 
menor média do Rn, o que pode estar relacionado 
com o início das queimadas na região, aumentando 
a quantidade de partículas na atmosfera (Biudes 
et al., 2009). 

Além do processo de absorção, ocorre na 
atmosfera o processo de difusão da radiação solar, 
determinada pelos seus diferentes constituintes. Esse 
processo é responsável por perdas apreciáveis da 
radiação que retorna para o espaço em virtude das 
dimensões das partículas determinantes da difusão. 
Essa difusão específica, determinada por partículas 
de maiores dimensões, não é dependente do com- 
primento de onda da radiação, havendo difusão 
para todos os comprimentos de onda e resultando 
a cor branca da luz difusa (Reichardt & Timm, 
2008). Esse fenômeno de nevoeiro, que é branco 
acinzentado, é muito comum no período seco. 

Os meses de fevereiro de 2007, janeiro de 


2008, janeiro de 2009 e novembro de 2009 apresen- 
tam linhas diagonais e estruturas de malha, carac- 
terísticas de sistemas com mais de duas frequências 
harmônicas, com períodos dentro de outros perío- 
dos. Nessas estruturas de malha (janelas de ciclos 
dentro de ciclos), em uma inspeção mais próxima, 
é possível observar estruturas de pequena escala, a 
sua textura, classificada em pontos isolados, linhas 
diagonais bem como linhas verticais e horizontais; 
ainda a combinação de linhas verticais e horizontais, 
obviamente formas retangulares de grupos de pontos 
de recorrência. Os descontínuos de pontos seguidos 
de retas paralelas à linha identidade explicitam o 
caráter oscilante entre aleatório e periódico dos 

r 

estados do sistema (Marwan et al.; 2007). E possível 
perceber faixas brancas características de dados 
não estacionários provenientes de estados menos 
frequentes, transições de fase no mecanismo do 
balanço de radiação. Observa-se ainda que o padrão 
que mais se destaca são as estruturas paralelas à 
linha identidade, o que se deve a forçante diária 
da energia solar. 

As Medidas de Quantificação de Recor- 
rência, que quantificam o número e duração das 
recorrências do sistema dinâmico representado 
pela sua trajetória no espaço de fase, oscilaram 
em estados transitórios sazonalmente (Tabela 2). 


Tabela 2 - Análise sazonal da taxa de recorrência (RR) e do de- 
terminismo (DET) em uma floresta sazonalmente inundada no 
norte do Pantanal. 



RR 

DET 

Seco 2006 

0,1476 

0,9521 

Chuvoso 

0,1412 

0,8985 

2006/07 



Seco 2007 

0,1713 

0,9916 

Chuvoso 

0,1463 

0,8824 

2007/08 



Seco 2008 

0,1699 

0,9807 

Chuvoso 

0,1486 

0,9163 

2008/09 



Seco 2009 

0,166 

0,9688 

Chuvoso 

0,1394 

0,8857 

2009/10 




A taxa de recorrência (RR) varia de 0 a 1 , 
indicando a probabilidade de que um estado retorne 
à vizinhança de raio s no espaço de fase (Webber 
Jr & Zbilut, 1994). Os resultados encontrados para 


Ciência e Natura, Santa Maria, v. 35 n. 1 jul. 2013, p. 086-097. 


95 



SAZONAL1DADE 


Figura 7 - Resultados sazonais da (a) taxa de recorrência e do (b) 


determinismo do saldo de radiação do Pantanal. 


RR (0,1439±0,0037 para os períodos chuvosos e 
0,1691 ±0,0022 para os períodos secos) indicam 
que, embora o sistema apresente uma caracterís- 
tica senoidal ligado à forçante da radiação solar, 
as trajetórias no espaço de fase não retornam 
exatamente à mesma vizinhança em seus ciclos 
diários (Figura 7). 

O determinismo DET é a razão entre o 
número de pontos pertencentes às diagonais e o 
número de pontos recorrentes. Logo, quanto mais e 
maiores forem as diagonais de um GR, mais pares 
de trechos do atrator se comportam de maneira 
semelhante, mostrando que uma previsibilidade 
maior é possível. Esses valores variam de 0 a 1; 
quando o valor tende a 0, o sistema é estocástico 
e quando tende alo sistema é determinístico 
(previsível) (Marwan et al., 2007). Os resultados 
encontrados para DET possibilitam afirmar que 
com um modelo apropriado é possível determi- 
nar (ou prever) o saldo de radiação no Pantanal 
com mais de 89% de certeza no período chuvoso 
(0,8957 ± 0,0133) e com mais de 97% de certeza 
no período seco (0,9733 ± 0,0147). Os valores de 
DET para o período chuvoso são limitados pela 
participação da nebulosidade na dinâmica do 
balanço de radiação; enquanto que muito possi- 
velmente para o período seco sofrem a influência 


dos aerossóis e emissões por queimadas (Artaxo 
et al.; 2005, Biudes et al., 2009). 

r 

E possível que os resultados obtidos tenham 
relação com o desmatamento no saldo de radiação 
das florestas sazonalmente alagadas do Pantanal. 
Contudo um número maior de observações e de 
séries temporais mais longas se fazem necessá- 
rias para se afirmar com segurança sobre esses 
possíveis efeitos. 

CONCLUSÃO 

As análises dos dados de saldo de radiação 
quanto às estimativas da dimensão de correlação 
e da dimensão embutida sugerem a existência de 
atratores estranhos de baixa dimensão fractal com 
sensível variação sazonal. 

O saldo de radiação no Pantanal segue uma 
dinâmica periódica acoplada a uma componente não 
linear complexa. Houve variabilidade sazonal do 
comportamento dinâmico entre os períodos secos e 
chuvosos das medidas complexas de taxa de recor- 
rência (RR) e de determinismo (DET). Observou-se 
que o ecossistema é bastante sensível às mudanças 
temporais pela presença de intermitências. 

Os GRs contribuíram para melhor caracte- 
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rizar a complexidade do sistema, exibindo informa- 
ções importantes sobre sua dinâmica e facilitando 
a interpretação desses estados em séries temporais. 
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